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Preface 


Surfing on the Ocean of Numbers - why this title? Because this little book does not 
attempt to give theorems and rigorous proofs in the theory of numbers. Instead it will 
attempt to throw light on some properties of numbers, nota bene integers, through a 
study of the behaviour of large numbers of integers in order to draw some reasonably 
certain conclusions or support already made conjectures. But no matter how far we 
extend our search or increase our samples in these studies we are in fact, in spite of 
more and more powerful technologies, merely skimming the surface of the immense 
sea of numbers. - Hence the title. 


Most books in Mathematics are used as reference books. I still consult my first 
Number Theory book which I bought in 1949 - Elementary Number Theory by 
Uspensky and Heaslet. It was when I was young and enthusiastic and dreamt about 
becoming a Mathematician. I am still enthusiastic but I became a Physicist instead.. 
However, [ stayed on the theoretical side and avoided to have to much to do with 
things that can break. But even so experiments are a major source of knowledge and 
maybe this book shows a little of a Physicist's approach to Mathematics. Most results 
are presented or supported by tables and graphs. All calculations have been carried out 
on a Pentium 100 Mhz laptop using Ubasic as a programming language. Finally, the 
author has tried to make a book which should be easy and pleasant to read. 


A word about the beauty of Mathematics and Number Theory in particular. The 
crystallized truth of a theorem, where a whole spectrum of mathematical thoughts 
come together to form an entity, is like a painting where designs and colours merge 
into a work of art. But sometimes it is not the finished result which is the most 
interesting - it could be the unsolved problem itseif. Why? Maybe it is the challenge of 
getting somewhere with it or the hours and days of thinking and trying that occupy the 
mind in a positive sense different from the problems of our time. It all brings piece to 
the mind - its like walking in the silence of the forest enjoying the trees, the sun and 
the blue sky , and should it happen that all the bits and pieces suddenly fall into place 
to give a solution then it is the most sublime experience for the human mind - eureka! 
But then the interest in the problem fades away unless solving the problem created 
new ones - and that is almost always the case. 


Most topics in this book have been selected from Only Problems, Not Solutions by Е. 
Smarandache. Others have bee suggested by Dr. R. Muller of Erhus University Press. 
Afew problems which the author has found. interesting originate from the Numbers 
Count Column of Personal Computer World. This journal has had great importance 
for the author as a source of recreational Mathematics and I take this opportunity to 
thank the Editor of this column Mike Mudge for all correspondence and 
encouragement he gave me in the past. 


Illustrations, graphics, layout and final editing up to camera ready form has been done 
by the author. Tables have been created by direct transfer from computer files 
established at the time of computation to the manuscript so as to avoid typing errors. 


This book has come into being thanks to R. Muller at Erhus University Press who has 
never failed an opportunity to give his support and encouragement. Rapid e-mail 
exchange between him in the USA and me in France has greatly facilitated our work. I 
also thank Dr. Muller's colleagues for their help. Many thanks are also due to my son 
Michael Ibstedt for his help and advice concerning computer equipment and software. 


Last but not least my warm thanks to my dear wife Anne-Marie for her encouragement 
and endless patience with a husband who does not always listen because his mind is 
Somewhere else. 


February 1997 
Henrv Ibstedt 
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Chapter I 


On Prime Numbers 


This chapter deals with some computational observations on prime numbers and their 
distribution. These computations only skim the surface of the ocean of integers but 
they give an idea of the general behavior of primes and often support some of the 
many conjectures that are made concerning primes. Computer programs have been 
written in Ubasic with extensive use of some of the built in functions of this language. 
In some cases the use of these functions will be illustrated with a few lines of program 
code. Most results are given in tabular and/or graphical form. 


1. On the Sequence a-pn + b 


In his book Only Problems, Not Solutions! Е. Smarandache asks the following 
question: 

If (a,b)=1, how many primes does the progression apn + b, where pn is prime and 
ne {1,2,...}, contain? 


Already for а=1 and b=2 we run into the classical unsolved problem "Аме there 
infinitely many twin primes? ". The answer to how many? is certainly equally difficult 
for other sets of parameters a,b. However, some interesting information on how a-pn + 
b behaves will be obtained for the first 10,000,000 primes pp. 


Let m be the number of primes produced by a-Pn + b for n<N, i.e. if ap; + b is prime 
we can write ар, + b = qm where qm is prime. The following Ubasic program lines 
have been used to determine whether а-рь + b is prime or not: 

while N«10000001 

p-nxtprm(p) 

inc N 

c=a%* pt+b% 

if nxtprm(c-1)=c then inc m 

wend 
The program has been implemented for a set of values of the parameters a and b. The 
result is shown in table 1. It is interesting to visualize the result. Because of the 
logarithmic behaviour of the distribution of primes it is reasonable represent m/N as a 


1 Unsolved problem number 17. 


function of logioN rather than as a function of М. For this reason the value of m has 
been recorded during the computation for N=10, 107, 103, 10* 10°, 109, and 107. 


Table 1. Number of primes m in the progression a-pr.-5 for n«N 





a,b /N 10 102 103 104 105 106 107 
1,2 5 25 174 1270 10250 86027 738597 
2,1 5 25 166 1221 9667 82236 711153 
3,2 8 47 290 2350 18919 160127 1392733 
4,1 3 2] 145 1108 9314 78676 685069 
Si 5 26 188 1492 12020 103010 903165 
3] 7 39 277 2175 18019 153925 1342255 
7.2 4 23 167 1288 10634 91232 


The graphs in figure 1 show m/N (y-axis) as a function of logicN (x-axis) where m is 
the number of primes of the form ap,+b for n<N. Figure 16 is an enlargement of figure 
la for large values of М, (N <10 ). The eight curves correspond to the following sets of 























Figure 1а. 
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pare 2 1215 listed in the order in which the curves appear from top to bottom in the 
right + 74 side of the two figures: (a,b)=(3,2), (6,1), (5,2), (7,2), (1,2) (2,1), (4,1) and 


Table 2. Number of primes m generated оу ратіб for NSN 

















b/N 10 100 i000 10000 700000 00000 
2 5 25 174 1270 10250 86027 
4 4 27 7 1264 10214 85834 
5 7 48 344 2538 20472 170910 
8 5 24 178 1308 10336 858% 
10 5 34 231 1682 13653 114394 
12 7 48 340 2515 | 20462 | 171618 
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Figure 2. The ratio m/M plotted against logN for b=2,4,6,8,10,and 12 
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The number of primes m in the sequence a-p,+b for n = 1, 2, ... N is illustrated in 
figure 2 fora = 1 and b = 2, 4, 6, 8, 10 and 12 where the ratio M=m/N is plotted 
against logN. The corresponding numerical results are given in table 2. 


The general appearance of the graphs for various values of the parameters (a,b) in 
ар:+Ы is very similar. In particular figure 1b shows an interesting picture of curves 
running parallel to one another and in particular to the one for (a,b)=(1,2), that is the 
curve for p+2 which corresponds to prime twins for which we have the classical 
conjecture that there are infinitely many. This makes the following conclusion 
reasonable. 


Conjecture: The progression ap«tb , (a,b)-1 contains infinitely many prime 
numbers. 


2. Prime Number Gaps 


Smarandache asked how many primes there is in the progression ари "р. For a=1 and 
b=2 the question is equivalent to ‘how many twin primes are there?'. Since we have a 
very stable conjecture that there are infinitely many we now want to know something 
about their distribution and also about the distribution of other prime number gaps 
£7pa«1-pn. With a small change in the Ubasic program used in the previous section we 
can study the distribution of primes over gaps g=2, 4, 6, ... 


p=3 

while p<N 

ЯР 

p=nxtprm(p) 

u=(p-qv2 

inc Ки) ‘Count the number of gaps = p-q. 

if p(u)-0 then рид “Store the smallest prime for which the 
wend “gap occurs in р). 


This program was run for primes р<№=2.10° The result is shown in table 3,where f is 
the number of gaps g and p the prime number for which the gap first occurs, N<2-10’. 
All gaps g<292 except 264, 278 and 286 are represented in the table which is 
arranged so that gaps g=2 (mod 6), g=4 (mod 6) and g=0 (mod 6) are found in 
separate columns . Gaps g=0 (mod 6) occur much more frequently than the other two. 
This is illustrated in figure 3 which also shows that In f as a function of g has a near 
linear behaviour for all three types. The “wild” behaviour fot gaps>250 would 
certainly correct itself if the range of primes in the study were extended. The area 
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below the curves for g«2 (mod 6) and g=4 (mod 6) are equal as will be shown shortlv. 
The curve for gz2 (mod 6) behaves уегу well while the one for р=4 (mod 6) shows an 
interesting ripple effect. In particular it shows а “bump” for g=70 which showed up 
already in the smallest sample N<10° for which g=70 first appeared. What causes this 
high frequency for g=70? 


For a prime number p25 we have pex! (mod 6). Let а and р be two consecutive 
primes forming the gap g=p-q. We distinguish between the following cases: 


Shift 
1. а=! (mod 6) and ре! (mod 6) = g=0 (mod 6) ++ 
2. q=! (mod 6) and р=-1 (mod 6) = g=4 (mod 6) +- 
3. 4=- (mod 6) and pei (mod 6) 2 2=2 (mod 6) -+ 
4. д=-1 (mod 6) and pe-1 (mod 6) => g=0 (mod 6) - 


А sequence of consecutive primes (with the first prime =5) can be characterized by the 
shifts: 
aH path tt HH tte, 


The longest sequence of consecutive primes = 1 (mod 6) for р<10? is of length 18: 


450988159, 450988177, 450983207, 450988231. 450988241, 450988261, 
45 50988297, 450988333, 450988339, 450988381, 450988387, 450988399, 
450988411. 450988423, 450988441, 450988471, 450988477, 450988567 


and the longest sequence of consecutive primes = -1 (mod ©) for p<10 is of length 22: 


19201, 766319231, 766319237, 766319249, 766319261, 
19291, 766319339, 766319357, 766319363, 766319369, 


755319189. 76 
7 


б 
766319273. 766 


663 


766319423, 766319441, 766319453, 766319483, 766319507, 766319549, 


ToC 19573, 765319279, 766319521, 756319027 





ities of finding the next shit to be +-, ++, -+ or — are 
НЕЕ, rf iti) and Ро бутан we would 
have F;-F470.25 and Ео=0.5. This is not the case. 
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Before looking into this let's first consider а related question: Do primes=I (mod б) 
(notation f.) occur with the same frequency as primes = -1 (mod 6) (notation f.) ? 
Table 4 shows a study of the number of primes f. and f. congruent to -1 respectively 1 
(mod 6) for primes less than 10* for k-1, 2, 3, ... 9. 


Within the range of this study we have Ё. > f +. However, the ratio г =(f . - К УЕ. - б.) 
is decreasing. Will eventually f . < f + ? 


We have proved that fz = f4 (assuming the last shift to be +-). We will now study the 
relative frequencies defined through 


Fo=fa/(fotfz+fa) Бе) —— Fo-fo/(fotfz*f4) 
Again we have F, = F, and of course Fo=1 - 2-Е». То study how F, varies as we 
increase the number of consecutive primes p«10* the execution of the program for gap 
statistics was stopped for k=1, 2, ... 9 to produce the data shown in table 5. 


Table 4. Number of primes = -1 respectively 1 (mod 6). г=(Е - №)/( + f+} 


End 
k 1i 2. 3 45 5 VRS 8 9 





f. 1 12 85 ` 616 4805 39264 332383 2880936 25424819 
1 


f+ 11 81 611 4785 39232 332194 2880517 25422713 
f.-f. 0 1 4 5 20 32 189 419 2106 
4f. 2 23 166 1227 9590 78496 664577 5761453 50847532 

ri О 435 241 4] 2] 4 3 0.7 0.4 
Table 5. Prime number gap distribution (mod 6) for primes «10* 

к  gexmod  ge4modé) ото lod Р; 
1 2 1 0 3 0.5 
2 9 8 7 24 0.354166667 
3 58 56 53 167 0.341317365 
4 379 378 471 1228 0.308224756 
5 2870 2868 3853 959] 0.299134605 
6 22839 22837 32821 78497 0.290941055 
7 189285 189284 286009 664578 0.284819088 
8 1616471 1616470 2528513 5761454 0.280566416 
9 14107250 14107249 22633034 50847533 0.277442162 
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2 f f f 

2 6388042 3 4 6386967 7 6 11407651 23 

8 5069051 89 10 6568071 139 8472823 199 
14 4690561 113 16 3527160 1831 6427670 523 
20 3528810 887 22 3030348 1129 4600962 1669 
26 2190452 2477 28 2386944 2971 4298663 4297 
32 1359889 5591 34 1430231 1327 2341569 9551 
38 1103677 30593 40 1308406 19333 1940894 16141 
44 796213 15683 46 687135 81463 1190342 28229 
50 678359 31907 52 511183 19609 856601 35617 
56 444581 82073 58 383239 44293 789454 43331 
62 247659 34061 64 253846 89689 466901 162143 
68 191321 134513 70 272834 173359 277514 31397 
74 137620 404597 76 122523 212701 233230 188029 
80 119756 542603 82 85030 265621 176328 461717 
86 63174 155921 88 65612 544279 133019 404851 
92 44723 927869 94 40821 1100977 71864 360653 
98 37946 604073 100 39504 396733 52752 1444309 


12 
18 
24 
30 
36 
42 
48 
54 
60 
66 
72 
78 
84 
90 
96 
102 
108 36484 2238823 
114 26413 492113 
120 23526 1895359 
126 14443 1671781 
132 8974 1357201 
138 6567 3826019 
144 4104 11981443 
150 4022 13626257 
156 2152 17983717 
162 1413 39175217 
164 661 20285099 166 553 83751171 168 1271 37305713 
170 607 27915737 172 430 38394127 174 729 52721113 
176 332 38069277 178 292 39389989 180 638 17051707 
182 238 36271601 184 235 79167733 186 342 147684137 
192 
198 
204 
210 
216 
222 
228 
234 
240 
246 
252 
258 
270 
276 
282 
288 


188 124 134065829 190 205 142414669 219 123454691 
194 109 166726367 196 112 7039633 221 46006769 
200 91 378043979 202 71 107534587 129 112098817 
206 44 232423823 208 56 192983851 141 20831323 
212 35 215949407 214 38 253878403 50 202551667 
In 21 327966101 220 36 47326693 31 122164747 
224 18 407866323 226 15 519653371 21 895858039 
230 17 07010093 232 3 525436489 23 189695659 
236 10 216668603 238 8 673919143 15 391995431 
242 8 367876529 244 5 693103639 7 555142061 
248 6 191912783 250 8 387096133 8 630045137 
254 3 1202442089 256 1 1872851947 2 1316355323 
260 3 944192807 262 } 1649328997 2 1391048047 
266 1 1438779821 268 1 1579306789 ү 649580171 
272 1 1851255191 274 1 1282463269 з 436273009 
284 2 1667186459 280 2 1855047163 2 1294268491 
290 1_ 1948819133 292 1 1453168141 
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100 


Figure 3. In(f) as a function of g for №2.107. се (mod 6) upper solid line, д=2 (mod 
6) lower solid line and gs4 (mod 6) dashed line. 


17 











Figure 4. Relative frequency of prime gaps . F2 = Fa (thick line) and Fo (thin line). 
Conclusion: Fs «0.45 and F = Е, 20.27 for primes <10°. 


The approach to the values 0.5 and 0.25 that one would have expected is very slow 
and is slowed down with increasing К, - one realizes that the interval 8<К<9 is ten 
times as large as the interval 0<k<8. The data used is cumulative but even if we 
consider only the interval between 10° and 10° we have F;7712,659,767/45,680,669— 
0.2771 compared to F7-0.2774 for the whole interval between 0 and 10°. 


Question: Given an arbitrarily small number 5>0, does a prime p; exist so that 
Е;<0.25+5 for all р»рі? 
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Chapter H 


On Smarandache Functions 


1. Smarandache - Fibonacci Triplets 
We recall the definition of the Smarandache Function S(n): 

S(n) - the smallest positive integer such that S(n)! is divisible by n. 
and the Fibonacci recurrence formula: 

Fa = Fai + Fa (n22) 
which for Fo = F; = 1 defines the Fibonacci series. 
We will concern ourselves with isolated occurrences of triplets n, n-1, n-2 for which 
S(n)-S(n-1 H S(n-2) and pose the questions: Are there infinitely many such triplets? Is 
there а method of finding such triplets which would indicate that there are in fact 
infinitely many of them? 
A straight forward search by applying the definition of the Smarandache Function to 
consecutive integers was used to identify the first eleven triplets [1] which are listed 
in table 1. As often in empirical number theory this merely scratched the surface of the 
ocean of integers. As can be seen from figure 1 the next triplet may occur for a value 
of n so large that a sequential search may be impractical and will not make us much 


wiser. 


Table 1. The first 11 Smarandache-Fibonacci triplets 


# n Sin Sirri Sin- 
1 11 11 5 2-3 
2 121 2-11 5 17 
3 4902 43 29 2.7 
4 26245 181 18 163 
5 32112 223 197 243 
6 64010 173 2.23 127 
7 368140 233 2.41 151 
8 415664 313 2.73 167 
9 2091206 269 2.101 67 
10 2519648 1109 2.101 907 
11 4573053 569 2-53 463 
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However, an interesting observation can be made from the triplets already found. 
Apart from n-26245 the Smarandache-Fibonacci Triplets have in common that one 
member is two times a prime number while the other two members are prime 
numbers. This observation leads to a method to search for Smarandache Fibonacci 
triplets in which the following two theorems play a róle: 


1. If n = ab with (a,b) = 1 and S(a) < S(b) then S(n) = S(b). 


П. If n = p° where p is a prime and a < p then S(p*)= р. 
































0 2 4 6 8 10 12 
Figure 1. The values of n for which the first 11 Smarandache-Fibonacci triplets occur. 


The search for Smarandache-Fibonacci triplets will be restricted to integers which 
meet the following requirements: 


n = хр? with asp and S(x)«ap (1) 
п-1 = уд? with bxq and S(y)«bq (2) 
n-2 = ді? with csr and S(z)«cr (3) 
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Table 2a. Smarandache-Fibonacci triplets 


N 


4 

11 

121 

4902 

32112 
64010 
368140 
415664 
2091206 
2519648 
4573053 
7783364 
79269727 
136193976 
321022289 
445810543 
559199345 
670994143 
836250239 
893950202 
937203749 
1041478032 
1148788154 
1305978672 
1834527185 
2390706171 
2502250627 
3969415464 
3970638169 
4652535626 
6079276799 
6493607750 
6964546435 
11329931930 
11695098243 
11777879792 


SIN) 
4 
И 


ж 


S(N-1) 


Table 2b. Smarandache-Fibonacci triplets 


N 


13429326313 
13849559620 
14298230970 
14988125477 
17560225226 
18704681856 
23283250475 
25184038673 
29795026777 
69481145903 
107456166733 
107722646054 
122311664350 
126460024832 
155205225351 
196209376292 
210621762776 
211939749997 
344645609138 
484400122414 
533671822944 
620317662021 
703403257356 
859525157632 
898606860813 
972733721905 
1185892343342 
1225392079121 
1294530625810 
1517767218627 
1905302845042 
2679220490034 
3043063820555 
6098616817142 
6505091986039 


13666465868293 


SIN) 
4778 


SIN-1) 
1597 
2474 
1847 
2986 
3118 
1823 

463 
1951 
8819 
3722 
6043 
6673 

10463 
2578 
4034 
3257 
1567 

11273 
2803 

12658 

18118 

20302 

10874 
3557 

13402 
10214 
12022 
9293 
5807 
8318 

21478 

7459 
12202 
20206 
19862 
16442 


p.q and г are primes. We then have S(n)-ap,S(n-1)-bq and $(п-2 )=ст. From this and 
by subtracting (2) from (1) and (3) from (2) we get 


ар = bq + cr (4) 
хр" - уд? = 1 (5) 
уф - zr =1 (6) 


The greatest common divisor (p*.q?)=1 obviously divides the right hand side of (5). 
This is the condition for (5) to have infinitely many solutions for each solution (p,q) to 
(4). Such solutions are found using Euclid's algorithm and can be written in the form: 
x = xo + qt, у= yo- рї (57) 
where t is an integer and (хо, vo) is the principal solution. 
Solutions to (5°) are substituted in (6^) in order to obtain integer solutions for 2. 
z= (уд? - 1yr* (6°) 
Solutions were generated for (a,b,c)=(2,1,1), (a;b,c)-(1,2,1) and (a,b,c)-(1,1,2) with 
the parameter t restricted to the interval -11<t<11. The result is shown in table 2. 


Since the correctness of these calculations are easilv verified from factorizations of 


S(n), S(n-1) and S(n-2) these are given in table 3 for two large solutions taken from an 
extension of table 2. 


Table 3. Factorization of two Smarandache-Fibonacci triplets 


n= 16,738,688,950,356=22-3-31-193.15,2692 5(п)= 2-15,269 

п-1 = 16,738,688,950,355=5-197-1,399-1,741-6,977 S{n-1)= 6977 

n-2- 16,738,688,950,354=2-72-19-23-53.313-23,561 5(п-2)= 23,561 

n= 19,448,047 ,080,036=22-32.432. 17,093? S(n)- 2-17,093 

п-1 = 19,448,047 ,080,035=5-7-19-37-61-761-17,027 5(п-1)= 17,027 

п-2 = 19,448,047 ,080,034=2-97-1,609-3,631:17,159 5(п-2)= 17,159 
що 

Conjecture: 


There are infinitely many triplets n, n-1, n-2 such that S(n)-S(n-1 y*S(n-2). 
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Questions: 


1. Itis interesting to note that there are only 7 cases in table 2 where S(n-2) is two 
times a prime number and that they all occur for relatively small values of n. 
Which is the next case? 

2. The solution for n-26245 stands out as a very interesting one. Is it a unique case 
or is it a member of a family of Smarandache-Fibonacci triplets different from 
those studied here? 


References: 


C. Ashbacher and M. Mudge, Personal Computer World, October 1995, page 302. 


2. Radu's Problem 


For a positive integer n, the Smarandache function S(n) is defined as the smallest 
positive integer such that S(n)! is divisible by п. Radu [1] noticed that for nearly all 
values of n up to 4800 there is always at least one prime number between S(n) and 
(8+1) including possibly S(n) and S(n*1) The exceptions are n-224 for which 
5(пу-8 and $(п+1)=10 and n-2057 for which S(n)-22 and S(n*1)-21. Radu 
conjectured that, except for a finite set of numbers, there exists at least one prime 
number between S(n) and S(n*1). The conjecture does not hold 1f there are infinitely 
many solutions to the following problem. 


Find consecutive integers п and п- 1 for which two consecutive primes рь and D-1 
exist so that рь < Міт(5(п),5(п- Ij) and рь-: > Махіб (т), 5(п- 1)). 


Consider 
nti = xp, (1) 


and < 
n= Ypi (2) 


where p; and розі are consecutive prime numbers. Subtract (2) from (1). 


хр? - урга = 1 (3) 
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The greatest common divisor (p;’,pr-1°)} = 1 divides the right hand side of (3) which is 
the condition for this diophantine equation to have infintelv many solutions. We are 
interested in positive integer solutions (x,v) such that the following conditions are 
met. 


5(п+1) = spr, i.e. S(x) < sp; (4) 

S(n) = spr: , Le. $(у) «sp: (5) 
In addition we require that the interval 

Spr < q < sp,-1° is prime free, i.e. that q is not a prime. 


Euclid's algorithm is used to obtain principal solutions (хо. vo) to (3). The general set 
of solutions to (3) is given bv 


X = Xo* pi^. y7yo-prt (7) 
with t an integer. 


These algorithms were implemented for different values of the parameters d-p.-i-p., 5 
andt. The result was a very large number of solutions. Table 4 shows the 20 smallest 
(in respect of n) solutions found. There is no indication that the set would be finite. 
One pair of primes may produce several solutions. 


Within the limits set bv the design of the program the largest prime difference for 
which a solution was found was 4-42 and the largest exponent which produced 
solutions was 5-4. Some numericallv large examples illustrating these facts are given 
m table 5. 


To see the relation between these large numbers and the corresponding values of the 
Smarandache function in table 5 the factorizations of these numbers are given below: 


1182293664715229578483018 = 2-3-89-193-431-16127812 
1182293664715229578483017 = 509-3253-16128232 

11157906497858100263738683634 = 2-7.372.5667 1.553335 
11157906497858100263738683635 = 3-5-11-192.16433.553373 
17549865213221162413502236227 = 3-112.307.12671-553333 

17549865213221 162413502236226 = 2-23-37-71-419-743-5533/3 
27032997592120525363470705182284857039 1314 = 2-33-47-1289-2017-119983- 1674414 
270329975921205253634707051822848570391313 = 37-23117.24517-38303-167443* 


Table 4. he 20 smallest solutions which occurred for 522 and d-2 








Ф n Sin) $i) р р» і 
1 265225 206 202 199 211 0 
2 843637 302 298 293 307 0 
3 6530355 122 118 113 127 -1 
4 24652435 926 922 919 929 0 
5 35558770 1046 1042 1039 1049 0 
6 40201975 142 146 139 149 1 
7 45388758 122 118 113 127 -4 
8 46297822 1142 1138 1129 1151 0 
9 67697937 214 218 211 223 0 
10 138852445 1646 1642 1637 1657 0 
11 157906534 1718 1714 1709 172] 0 
12 171531580 1766 1762 1759 1777 0 
13 299441785 2126 2122 2113 2129 0 
4 551787925 2606 2602 2593 2609 0 
15 1223918824 3398 3394 3391 3407 0 
16 1276553470 3446 3442 3433 3449 0 
17 1655870629 3758 3754 3739 3761 0 
18 1853717287 3902 3898 3889 3907 0 
19 1994004499 3998 3994 3989 4001 0 

20 2256222280 4166 4162 4159 4177 0 
Tabie 5. Four numerically large solutions 
Pairs of consecutive integers S) d s 1 Dy pet 
1182293664715229578483018 3225562 42 2 -2 1612781 
1182293664715229578483017 3225646 1612823 
11157906497858100263738683634 165999 4 3 0 55333 
11157906497858100263738683635 166011 55337 
17549865213221162413502236227 165999 4 3 -i 55333 
17549865213221162413502236226 166011 55337 
270329975921205253634707051822848570391314 669764 2 4 O 167441 
270329975921205253634707051822848570391313 669772 167443 





It is also interesting to see which are the nearest smaller Р; and nearest bigger P,-1 
primes to S; = Min(S(n),S(n+1)) and $; = Max(S(n)S(n*1)) respectively. This is 
shown in table 6 for the above examples. 


Table 6. Pk < $1 < $2 < Рея 





Рк З S2 Рен С = Pu - Рк 
3225539 3225562 3225646 3225647 108 
165983 165999 166011 166013 30 
669763 669764 669772 669787 24 
Conclusion: 


There are infinitely many intervals {Min(S(n),S(n+1)),Max(S(n),S(n+1))} which are 
prime free. 
References: 


LM. Radu, Mathematical Spectrum, Sheffield University, UK, Vol. 27, No. 2, 1994/5, 
p. 43. 


3. The Smarandache Ceil Function 


Definition: For a positive integer n the Smarandache ceil function of order k is 
defined through ' 


Sx(n) = m where m is the smallest positive integer for which n divides m“. 


In the study of this function we will make frequent use of the ceil function defined as 
follows: 


[х1 = the smallest integer not less than x. 


The following properties follow directly from the above definitions: 


1. Si(n)=n 
2. Sp") = р for any prime number р. 


3. For distinct primes р, 4, ...г we have S(p*q^ 55) = pog PM Lu Sei 


Theorem L Sx(n) is a multiplicative function. 


! This function has a great resemblance to the Smarandache function. It’s definition was proposed by 
K. Kashihara (Japan) and conveyed to the author by R. Muller. 
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A function Кп) is said to be multiplicative if for (плз) = 1 if it is true that f(nim) = 
f(ni)f(n;). In our case it follows directly from (3) that if (nı,n2) = 1 then Зп») = 
Si(n1)Sk(n2). 


However, consider n-nim when (nı,n2)¥1. In a simple case let n=m,p* and пусту p? 
with (m,,m2)=1 we then have Si(n) = Simi)Sm;)p ©" which differs from 
Si(n1)Si(n2) whenever [to* yk | + [ak] + [Вк]. In fact one easily proves that 
[(o*BYkl- Го + Грк Т ог [(а+В УК] = Галю 1 + ГАК - 1. 

Theorem П. S,.;(n) divides Sy(n) 


Express п in prime factor form п=р°д?.-т? and apply (3). We then see that all prime 
powers in Sy.,(n) are less than or equal to those of Зк(п), i.e. Sx-1(n) | Зк(п). 


Theorem Ш. For sufficiently large values of k we have Ѕұ(п)= Пр where the product 
is taken over all distinct primes pi of n. 


By extending the argument in theorem II we have that, if j=max(a, В, ... 8) then 
Si(n)-pq--1 for k 2 j. 


Corollary 1. Sp) = p for any prime number р. 

Corollary 2. If n is square free then S2(n) = n. 

Theorem IV. К exists so that и (п!) = рй, where р is the largest prime dividing п. 
p? denotes the product of all primes less than or equal to p. Let's write n! in prime 
factor form. n!=2°3°..-p’, where a>B>....>y. In order to apply theorem Ш we need to 
find a. Consider 1-2-3-4-5-6 ... n. This product contains [п!/2] even integers, [п!/4] 


multiples of 4, etc ....and finally [n'/2^] multiples of 2°, where 25«n'«2*. $ is 
determined by ó-[log n/log 2]. From this we find that S,(n!) = p# for 
б 


к=а = Y [n!/2"] 


r-1 
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Table 7. The Smarandache сей function 


оо кю ООо ьо рм 


З З о o 


15 
12 
7 
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Calculations. Calculation of S,(n) for п < 1000 were carried out in Ubasic, which has 
a built in ceil function. The result is shown in table 7. Since Sn) = n for square free 
numbers these have been excluded from the table. When S;(n) is square free the 
entries for larger values of k become repetitive. Instead of repeating these values the 
corresponding spaces in the table have been left blank. 


4. The Smarandache Pseudo Function Z(n) 


Definition’: Z(n) is the smallest positive integer m such that 1+2+...+m is divisible 
bvn 


Alternative formulation: For a given positive integer n, Z(n) equals the smallest 
positive integer m such that m{m+1)/2n is an integer. 


The following properties follow directly from the definition: 


. А 

. 203 

For апу odd prime number p, Z(p)-p-1 

. Ву extension of (3) we have Z(p5-p*-1 

. In the special case n=2* we have Z(25)-2* 1-1 


ль м м 


Calculation of Z(n) 


We need to find m so that m(m-*1)-2nk has a positive integer solution for the smallest 
possible positive value of k. 


-L-414-8kn 
m = 739 


For a given value of n the smallest square 1+8kn is found by executing the following 


program lines in Ubasic where effective use of the ISORT(x) has been made: 


2 Definition by K Kashihara (Japan) conveyed to the author by В. Muller, Erhus University Press, 
USA. 


30inck 

40 х=1+8*К*п 

50 if x<>(isqrt(x) 2 then goto 30 
etc - to evaluate m 


The complete program has been implemented for n<1000. The result is displayed in 
table 8. 


Theorem: If n=pq, where p and q are two distinct primes with g=q-p, then 
Z(n)}-Min(p(qk+1 Vg where pk+1=0 (mod g), q(pk-1Yg where рк-1=0 (mod g)) 


Proof: 

We have to consider three cases: 

1. p | m and а | (m*1) which, since we assume q>p, we distinguish from 
2. pl(m+1) and qlm 

3. pql(m+1) 


Case 1. Consider px=m and qy=m+1 which together with g=q-p gives 
р(х-у)=ру-1 (1) 


Since we must have р | (gy-1) we can put gy-1=pk where К | (х-у). Our solution for y 
then becomes 


y=(pk+1)/g with pk+1=0 mod g) (2) 
Inserting this in (1) results in 

P(x pk+1 Увуєрк 
from which 

x-(qk*1Yg (3) 


which we insert in т=рх to obtain 


m-p(qk-*1Yg where К is determined through pk+1=0 (mod g) 
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Table 8a. Z(n) for n < 1000, n non-prime 


— 0 


7 
3 
5 
8 
4 
8 
7 
5 
] 
8 
5 
6 
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Table 8b . Z(n) for n < 1000, n non-prime 
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Table 8c. (п) for n < 1000, n non-prime 





Table 8d. Z(n) for n < 1000, n non-prime 
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Case 2. Consider рх=т+1 and qy=m together with g=q-p. Repeating the calculation in 
case 1 results in: 


m-q(pk-1Yg where К is determined through pk-1=0 (mod g) 

Case 3. Consider m-pq-l. We must have pq-l«p(qk*lYg to beat case 1. This 
inequality can be written p(q(g-k-1)«g where g>k and q>2 from which we see that 
this inequality is impossible. The argument is similar to reject case 3 in comparison 
with case 2. m is therefore given by the minimum value for m as calculated in cases 1 
and 2, i.e. 

Z(n)-Min(p(qk*1Yg where pk+1=0 (mod р), q(pk-1Yg where pk-1=0 (mod g)) 
Corollary: 7(п) is not a multiplicative function. 
The above theory has been used to calculate Z(pq) for a few prime number gaps g=q- 


p. The result is shown in table 9. 


Table 9a. 2(п) for the first 10 prime gaps g=10 and g=30 


20711 | 18593119 8056874 
34571 | 23483491 782621 
60491 6024 | 28793731 10557522 
82919 24904 | 33223471 12182131 


116939 35047 | 42327811 
171371 17178 | 48052399 
181451 18102 § 52823599 
304679 91348 | 60434851 
338699 101551 | 63648259 
404471 40383 | 70626991 
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Table 9b. Z(n) for the first 10 prime gaps g=20 and g=40 


804509 
2712509 
9603701 
11716829 
15657749 
32798429 
34963469 
35963909 
46416869 
58660181 


949460) 
4321510} 
1757695] 


26396698 


37 


19373 
20849 
22613 
25301 
33287 
38501 
45053 
48947 
52361 
60209 


374538209 
433846841 
510443249 
639128561 
1106692889 
1480786961 
2027970689 
2393850929 
2739579881 
3622715321 


28090849 
97615018 
38283808 





Chapter Ш 


Loops and Invariants 


In many practical as well as theoretical processes we repeat the same operation on an 
object again and again in order to arrive at a final result, or sustain a certain state or 
maybe simply to see what is going to happen. Each time a hammer hits a nail the nails 
sinks a bit deeper until it can sink no further. This repeated operation has resulted in 
an invariant state. A different situation occurs in an engine when energy is used to 
make a piston perform cycles or loops. Many similar phenomena occur when repeating 
processes on numbers through iterations. Before going into the problems of this 
section let us consider the iteration process itself. 


Let (п) define an operation to be carried ош on n. If we apply this operation to I(n) 
itself we say that we perform an iteration and could write this as К(п)). After a 
number of iterations we could have something like I(I(...I(n)...)), which we will write 
I(x(n) to indicate the kth iteration. Alternatively we can use nx to denote the result of 
the kh iteration and noto denote the starting value. We then have 


nı = I(no), m = Ini), Фм? пк = Кик) 


Let us apply this to a simple case where I(n) is defined through Кп) = 1/n. If we take 
п0=1 then ny-1 for all К, i.e. the result of the iteration is invariant. If we take ш=2 
then n.-1/2 for odd values of k and пк=2 for even values of К, i.e. we have an iteration 
resulting in a loop of length 2. If we apply the iteration process to I(n)=n’ then the 
result will be forever increasing and we say that the iteration is divergent. 


Since we are dealing with a very important concept which has attracted a lot attention 
some of the topics in this section have been dealt with before. In particular Perfect 
Digital Invariants [1], which has recently been reactivated under another name 
"Steinhaus' problem" [2]. However, all results presented here have been generated in 
recent studies carried out by the author and have been retained even though some of 
them may duplicate earlier results. This proved necessary in order to arrive at a 
consistent presentation and some new results. J. S. Madachy has made me aware of 
more literature on this interesting topic some of which is listed in the references. 
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1. Perfect Digital Invariants and Related Loops 
For an arbitrary positive integer 
№ = ал10® + an110™ +... a110 + ао 
we define №; = (М) through 
(№) = an? + ал19 +... + а + ao", where q e N, 422. 


For a given positive integer по the iteration process n; = К), п; = Қт), ... is 
continued until one of the following situations 1s reached: 


І I(nx) = nx Perfect Digital Invariant (PDI) 
П. о) = nx Loop of length L 
Ш. (n) > х Divergent 


Case III is impossible. It. follows from the following determination of the upper search 
limit for each value of q. Let us consider the largest possible n-digit number М = 
999...99 = 107-1. We have IN) = n-9* and need to determine Nmax so that N-I(N)>0 
for all N>Nmax. Let u be the largest value of n for which 10"-1-n-9? < 0, i.e. 


10-1-u9*«0 while 10°1-1-(и+1)9%>0 
then there is a smallest positive integer a, 1<а<9, such that 
а-10° + 10" -1-а9 - ц-99> 0 


This gives Nmax = (а+1)-10* -1 as an upper limit for solutions. Nmax could be improved 
bv looking for a smaller value than 9 for the second most significant figure but this 
would give more complications than benefits in computer implementation. That Nmax 
exists proves that the iteration process does not diverge since after a number of 
iterations larger than Nmex a previously assumed value must be repeated completing a 
loop or collapsing on an invariant. 


Only а small subset of all integers < Nmax needs to be used as input numbers in a 


search program. The following two input criteria greatly reduce the computer 
execution time. 
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1. The order in which digits occur in an input number is of no importance. № = 
2337 will give the same result as No = 3732. A number whose digits are a 
permutation of an already used input number will therefore be rejected. As an 
example take q=5 for which Nmax = 299999. In this case an input number needs 
to have maximum 6 digits of which there can be at most 6 ones or twos and a 
maximum of 5 of the other digits. 


2. Input criterion number 1 is used so that the search always proceeds from a 
smaller to a larger input number. If for any input number No we have I(No)<No 
then the iteration process is aborted since I(No) has either been dealt with before 
or doesn’t meet criterion number 1. 


Complete solutions were calculated for q = 2,3,4, ... 15. Apart from the trivial case 
I(N)=1 these solutions are given in table 1 together with the upper search limit for 
each q. The longest loop is of length 381 for q = 14. There are no Perfect Digital 
Invariants (PDIs), i.e. solutions to N = I(N) for q=2, q=12 and q=15. 


References: 
1. Lionel E. Deimel Jr. and Michael T. Jones. Journal of Recreational Mathematics, pgs 
87-108, Vol. 14 .2 


2. Personal Computer World, page 333, January 1996 
3. Dean Morrow, Journal of Recreational Mathematics, pgs 9-12, Vol. 27.1 


40 


Table 1a. PDis and related loops 


q Nmax Length of Smallest Largest N 
S-EN Loop tem term 

2 199/4 8 4 145 3 

3 2999 ] 153 17 

1 370 16 

$=76 1 371 19 

1 407 2 

2 136 244 ] 

2 919 1459 11 

3 55 250 6 

3 160 352 3 

4 29999 1 1634 1 

1 8208 8 

S-153 1 9474 2 

2 2178 6514 6 

7 1138 13139 135 

5 299999 1 4150 1 

1 4151 1 

$=345 1 54748 10 

1 92727 1 

1 93084 1 

1 194979 1 

2 58618 76438 23 

2 89883 157596 1 

4 10933 73318 9 

6 8299 150898 13 

10 8294 183635 24 

10 9044 133682 33 

12 24584 180515 °з 

22 9045 167916 112 

28 244 213040 21 

6 3999999 1 548834 2 

2 63804 313625 5 

$=401 3 282595 845130 71 

4 93531 650550 5 

10 239459 1083396 167 

30 17148 1758629 150 
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7 


8 


9 


Nmax 
$=# М 
4107-1 


$=1012 


4.108-1 


5=1544 


4.109-1 


5=5058 


Table 1b. PDIs and related loops 


Length of 
Loop 


QNNM 


154 


Smallest 
term 


1741725 
4210818 
9800817 
9926315 
14459929 
2755907 
8139850 
2767918 
2191663 
1152428 
922428 
982108 
253074 
376762 
86874 
80441 
24678050 
24678051 
88593477 
7973187 
8616804 
6822 
146511208 
472335975 
534494836 
912985153 
144839908 
277668893 
180975193 
409589079 
52666768 
20700388 
62986925 
180450907 
42569390 
52687474 
322219 
41179919 
32972646 
2274831 


6586433 
9057586 
8807272 
9646378 
14349038 
16417266 
14600170 
18575979 
19210003 
19134192 
14628971 


77124902 
149277123 
153362052 


1043820406 
756738746 
951385123 

1339048071 
574062013 

1212975109 
931168247 
857521513 

1001797253 
708260569 

1298734342 

1202877221 

1724515947 

143071763] 
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Table 1c. PDIs and related loops 


Length of 
ioop 
1 


42 


Smallest 
term 
4679307774 
304162700 
123148627 
1139785743 
62681428 
20818070 
192215803 
32164049650 
32164049651 
40028394225 
42678290603 
44708635679 
49388550606 
82693916578 
94204591914 
4370652168 
2491335968 
4517134494 
6666140097 
416528075 
2181207047 
9005758176 
3967417642 
12650989279 
2075164239 
195493746 
101858747 
739062760 
872080538 
8922100 
98840282759 
2700216437 
4876998775 
16068940818 
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Largest 
term 


344050075 
7540618502 
9131926726 

13957953853 
15434111703 
14230723551 


11346057072 
71768229638 
33424168842 
36704410767 
103153306403 
28167146357 
71727610926 
98110415227 
128870085703 
127554589656 
106744983639 
134844138593 
169812860326 
165906857819 
176062673167 
785119716404 
1157645923834 
1281243062759 
12006154801 66 


100 
486 
1205 
1075 
1015 
1566 
609 
1142 
128 
1557 
1883 
5897 


Table id. PDIs and related loops 


а Nmax Length of Smallest Largest N 
525 N Loop tem term 

13 4.1013-1 1 564240140138 10 

3 3396705890823 6294418483143 240 

$=35177 5 2310198454262 10830899644190 50 

5 3656948275943 6405584099531 1601 

8 688271525816 5647840775906 8846 

9 1113928853354 5840462013812 503 

13 672769417360 6842950405261 240 

16 662218816395 10829615029038 3403 

18 100555910945 7735617503306 404 

22 221192011997 7722486519974 821 

31 759039131807 12916172647637 625 

100 98183216073 10719727658703 6794 

146 111240729637 11270705761846 11753 

14 4.1014-1 1 28116440335967 60 

4 6674458190799 74548238736415 9 

$=40180 4 27510477911590 47800729611562 428 

5 841332967215 57268678554888 154 

6 5911230616470 5591860436027 1 441 

6 22955961974580 94220062144011 311 

7 5833130055708 97267770775241 35 

14 1447300158177 92876091 448633 $119 

15 4423275478678 78789825354783 1698 

65 931899363208 118800319349172 3517 

96 157270047611 92354596687594 9448 

381 78924999008 142415666495594 18959 

15 4.1015-1 3 255349823145519 447090882837 630 17 

8 76058866219899 729415146625008 10157 

$=48785 11 70433029388274 657638038056753 518 

12 41845353013296 632509424234352 1261 

15 14745577655668 619115571288208 177 

19 14746637048586 895875047683584 362 

30 1096193499692 1030397118565900 941 

46 10059036985836 662536011243090 1302 

75 593343453343 864468451964287 18270 

216 92683904599 1065173414568230 2400 

362 595476555320 1145004602355100 13379 
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2. Тһе Squambling Function 


The squambling function U(N) originates from The Penguin Dictionary of Curious 
and Interesting Numbers, page 169. In its present form the problem was formulated by 
С. Samson in the November 1995 issue ої Personal Computer World: 


For a given integer N = ао + a110 +... tan10" > 10", an = 0, the squambling function 15 
given bv 


U(N) = a? + ат! +... + as. 


Iterating the squambling function will result in an invariant or a loop. As in the 
previous problem there are no starting values for which the process diverges. To see 
this consider 


U(N) = ат + a e ар «91 99+... 9*7? = 819" -1y8 


If № exists for which U(N) < № for all N > № then the squambling function must 
return to one of its previous values after a sufficiently large number of iterations, ie. 
the process converges (or goes into a loop). No exists and is estimated to be less than 
10“ from the inequality below 


U(N) < 81(977-1y8 < 10° < N 


More detailed considerations mav considerably lower this upper limit for No, but in 
any case the upper limit is so large that a complete search for invariants and loops is 
hardly feasible. The result of a search for squam-invariants and squam-loops for 
integers up to 10° is summarized in table 2, where L=the total number of numbers in 
the loop (-length of the loop, which equals 1 for an invariant), M-the smallest 
number in the loop, N-the largest number in the loop and Q-the number of loops 
found for initiating integers less than 10*. 


Table 2. Squam-loops and squam-invariants 


2077 
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If there are more invariants and loops then they must have a smallest element greater 
than 105, or more precisely a smallest element which is larger than the upper limit for 
a previous search. This has been used in an extended search for starting integers up to 
50,000,000. No new loops were found. Question: Are there any? 


Let us finally consider what happens if we reduce the powers to which each digit is 
raised by one. This makes the process a fortiori convergent as can be seen from the 
arguments used previously. The iteration process was examined for starting integers 
up to 10°. No loops were found. Apart from iterations ending on one-digit invariants a 
number of other invariants were found. 

Table 3. The number of start integers © below 105 resulting in an invariant I 


га | вә | 135 | 175 | 58 | 598 | 1306 | 1676 | 2427 | 
де | 347 | 492 | 54 | 319 | 128 | 102 | 20 | 27 | 


Question: Are there no loops in this case? 


3. Wondrous Numbers 


This study deals with the extended Wondrous Numbers Conjecture stated as follows 
by B.C. Wiggin [1]: 


Any integer п > (D-1) may be directed through a series of iterations as a 
function n = В (mod D). Starting with по the series of iterated Wondrous 
Numbers aren), N2, ... Пк, ... 
Definition: Let пк=Вх (mod D) 

if Ry-0 then n. =n D 

if ISRxSD-2 then ny-;7n(D- 1)-Rk 

if Ry-D-1 then ny-1=m(D+1 +1 


It is conjectured that the series ultimately converges to n«D. 


This study will show that the above conjecture does not hold. 


To examine the behaviour of the above process assume that for a given nx we have 


пу = 0 (mod D) 


Пк! = Rei (mod D) 
We have to consider three cases 


Case I 


Ку. =0 


Пк-2 = nk /D; п. = n/D? 
Case П 
0 < Re < D-2 
then 
Лео = Dk (D+1) - Rea 
From this we see that пк+ = 0 (mod D) and substituting from (1) that 
ngo < 01 *1/D) 
If Вазон stays within the above limits for І- 1, 2, ... m then 


Diam < m(1 + 1/0)" 


Case III 
Res ар 


As in case П m» =0 (mod D) but 
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п: < (1 +1/D) + 1 (4) 


For the series to diverge the case R=0 must occur with a frequency which is much 
higher than the expected 1/D (on the average). In a worst case scenario let us assume 
that R ¥ 0 so many times that it balances the effect of the occurrence of case I. 
Ignoring (4) which can be considered as compensated by ignoring Ку. in (3) we then 
have. 


ni + DYD = пі 


For D = 11 this gives m = 56 just to hold the balance. If remainders К are evenly 
distributed then m would varv around 10. It is very safe to say that the series does not 
diverge. However, it does not always converge to n « D as conjectured. Brendan 
Woods [2] reported that she failed to get termination with (Д.М) = (13,70), (14,75), 
(58,59), (82,83) and (198,199). The reason is that these cases produce loops as shown 
in table 4. № = 199 , D = 198 produces а very long loop of length 2279 requiring 2499 
iterations. 


Charles Ashbacher [3] investigated all divisors in the interval 3 < D < 12 for all initial 
values no < 14000000. Even with the most effective programming and up to date 
equipment this is a very impressive piece of work. In all cases iterations terminated on 
a one-digit number. He listed those cases for D-11 and D=12 where more than 1000 
iterations were required before the terminal value was reached and made the 
following observation "..for certain series of ascending n the number of cycles 
descended in steps of 2 contrary to the expected behaviour that the number of cycles 
rises with the size of the number”. He adds that no justification has been found for this 
and challenges readers to further explore this behaviour. The remainder of this section 
will be devoted to an explanation of this mystery though an analysis of the case D = 
11. A similar analvsis to the one below has been carried out for D = 12 with similar 
results. 


Since all iterations result in а one-digit terminal value all cases with по < 14000000 
which require more than 1000 iterations can be classified according to their terminal 
value. This is done in the column marked p in table 6, which contains Ashbacher’s 
table as a subset. There are only three different terminal values 3, 7 and 8. 


Table 5 shows the 180 first iterations for no = 1345680, which is the first entry in 
Ashbacher's table. Table 5 confirms the above theory that at least every second R is 
Zero. This explains the step 2 decrease with ascending n since this causes the terms to 
oscillate. In general every second term was out of bounds of the investigation which 
was limited to ng«14000000. These cases with long cycles are due to long “freak” 
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oscillate. In general every second term was out of bounds of the investigation which 
was limited to no<14000000. These cases with long cycles are due to long “freak” 
sequences if iterations with remainders R=0. Let's apply (3°) to the first occurrence of 
R=0 for even К in table 5. n; = 16148154, nig; = 167493499. 


Піст < 16148154(1 + 1/11)/121 = 167494096 
This, as is seen, is a very good approximation. 


Table 4. Loops instead of anticipated convergence 


n no Loop starts First term of Loop Length of 
і for: і finishes for: і 
13 70 92 1911 162 71 
14 75 13 166 91 79 
58 59 19 4002 555 537 
82 83 153 13120 925 773 
198 199 221 61380 2499 2279 


We also see from table 5 that it contains a sequence of numbers which also occur in 
table 6, i.e. the initiating value no = 1345680 is the “grandparent” of a whole family of 
lerger initiating values having the same terminal value 8. However, table 6 also 
contains initiating values with terminal value 8 which do not occur in table 5. At some 
stage, however, these will merge with the iteration process for the "grandparent" . 
These cases have been identified and labeled in column c (c-"child") in table 6. These 
"children" have been used as starting numbers for iteration processes to see at which 
point they will join the “grandparent” iteration. The result is shown in tables 7 and 8, 
where K-c is the number of iterations for the child (which may have a number of 
children of its own as can be seen from table 6), and K-p is the number of iterations 
for the grandparent before merging occurs. It is amazing how soon this happens. Only 
the 8-family child 11700624 makes it on its own almost to the end. N-max is the 
maximum value for the iteration process which occurs after K-m iterations. The 3- 
family is childless. 


References: 
l. B.C. Wiggin, Journal of Recreational Mathematics, pgs 52-56, Vol. 20.1 


2. M. Mudge, Personal Computer World, page 335, Dec. 1995 
3 С. Ashbacher, Journal of Recreational Mathematics, pgs. 12-15, Vol. 24.1 
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Table 5. The first 180 iterations for 0211, N=1345680 


ж 
»m 
3 
3 
x 
я 
5 
x 
җ 


16148154 
17616159 
19217627 
20964680 
22870551 
24949683 
27217828 
29692168 
32391447 
35336125 
38548499 
42052901 
45875885 
50046414 
54596080 
59559357 
64973843 
70880557 
77324236 
84353709 
92022227 
100387881 
109514053 
119469867 
130330761 
142179004 
155104367 
169204761 
184587007 
201367639 
219673784 
239644119 
261429949 
285196307 
311123241 
339407167 
370262365 
403922574 
440642807 
480701243 
524401350 
572074195 
624080941 
680815564 
742707889 


810226780 
883883759 
964236823 
1051894712 
1147521496 
1251841624 
1365645402 
1489794977 
1625230882 
1772979142 
1934159062 
2109991697 
2301809125 
2511064501 1 
2739343090 
2988374279 
3260044667 
3556412365 
3879722572 
4232424625 
4617190501 
5036935090 
5494838272 
5994369018 
6539311647 
7133794525 
7782321294 
8489805049 
9261605507 
10103569641 
11022075966 
12024082873 
13117181308 
14309652335 1 
15610529814 1 
17029668887 1 
18577820602 
20266713379 
167493499 
182720180 
199331105 
217452114 
237220489 
258785989 
282311986 
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0 
0 
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Table 6. Some Wondrous Numbers and their Genetic Relations, О=11 


1345680 
1468014 
1601469 
1747057 
1905880 
2038440 
2079141 
2223752 
2268153 
2425911 
2474348 
2597542 
2646448 
2699288 
2833682 
2887034 
2944677 
3091289 
3149491 
3212375 
3372315 
3435808 
3504409 
3606876 


4170535 
4292479 
4378181 
4460612 
4549674 
4682704 
4776197 
4866122 
4963280 
5108404 
5210396 
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7892741 
8050325 
8050330 
8201894 
8365657 
8610263 
8782172 
8782176 
8782178 
8947520 
9126171 
9393014 
9580551 
9580555 
9580557 
9760931 
9955823 
10057558 
10246920 
10246924 
10451510 
10451514 
10451516 
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1148 
1095 
1203 
1150 
1146 
1093 
1201 
1148 
1144 
1091 
1199 
1146 
1142 
1089 
1197 
1144 
1140 
1087 
1195 
1142 
1142 
1138 
1085 
1193 
1140 
1140 
1140 
1136 
1083 
1191 
1138 
1138 
1138 
1134 
1081 
1193 
1189 
1189 
1136 
1136 
1136 








10648288 
10860897 
10971881 
11178458 
11178462 
11311510 
11401647 
11401651 
11401653 
11616314 
11700624 
11735002 
11735026 
11848251 
11969324 
12194681 
12194685 
12339829 
12438160 
12438164 
12438167 
12438192 
12672342 
12764317 
12801820 
12801846 
12925364 
13057444 
3303288 
13303292 
13461631 
13568901 
13568906 
13568909 
13568928 
13568936 
13824373 
13924709 
13965621 
13965648 
13965650 
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Table 7. The 7-family. Grandparent = 2038440 
т. 


# — N-chid N-childone М-рсгепі K K N- N-max K 
one step  onestep  - - merger 
before before c p m 
merging merging 
1 3606876 293535231 2223752 46 3 26685021 8601368512990 562 
2 1024692 381088411 2887034 28 9 34644401 8601368512990 538 
3 1173500 764423308 5791086 97 25 69493028 8601368512990 591 
4 1173502 293535231 2223752 75 3 26685021 860136851 2990 591 
5 1396564 293535231 2223752 7 3 26685021 8601368512990 587 


Table 8. The 8-family. Grandparent = 1345680 
———Ш— MM EE E E ERE 
K 


# N-child  N-childone N- K N- N-max K 
one step опе - - merger - 
before before c р m 
merging merging й 
1 259754 504634735 3822991 66 25 45875885 22512799837 489 
2 805033 326613848 2474348 30 15 29692168 22512799837 463 
3 878217 388697375 2944677 32 19 35336125 22512799837 461 
4 100575 504634735 3822991 91 25 45875885 22512799837 514 
5 117006 38236 290 1048 1098 3476 43996530071824 348 
6 124381 177629694 1345680 6 1 16148154 22512799837 453 
7 135689 211393897 1601469 8 5 19217627 22512799837 451 


4. Iterating d(n) and o(n) - Two problems proposed by Е. Smarandache 


(a) Let n be a positive integer and (п) the number of positive divisors of n including 
1 and n. Find the к k for which d(d(... d(n)...)) 2 after k iterations, i.e. 
find k so that da(n) = 


d(n) is an important arithmetic function. We will look at its most important properties. 
The factors of p^, where p is a ише, are 1, p, p^, ... p*. Consequently Фр") = на. 
The number of factors in n = p^pP is easily seen to be (1+aX1+8) from which the 
following important theorem follows: 





! In fact k is a (single-valued) function of n. 
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Ш nı and n; have по common divisor, ie. are relatively prime which we write 
(ninj-1l, then d(ninj)sd(ni)d(n)) We say that the arithmetic function аф) is 
multiplicative. 
With n written in standard form n=p:°p2"...p," we have 

(п) = (1+0Х1+р)...(1+т) (1) 
We can now state: 

d(n)«n for all n. 

d(n)=1 if and only if n=1. 

d(n)=2 if and only if n is a prime number. 
From the above properties we see that d(n) is a measure of how far n is from being a 
prime number. The larger the number of factors of n the larger is (п), d(n) being 
equal to 2 only when n is a prime. This makes it interesting to try to answer 
Smarandache's question [1]: How many iterations k are required т daj(n)=2? 


Before looking at this problem let's make an important observation: 


Given an arbitrarily large positive integer k we can always construct 
infinitely many integers n for which dy)(n)>2 for all i<k and daj(n)=2. 


Here is how: Let pi, р», ... px be odd primes (not necessarily distinct). Make the 
following series of constructions: 


n, = р d(n;, =2 
пері" d(nj = pi 

n = pF "=p"  d(m)= р 
n, = pe d(n,) =n, 
n, = pus d(n,) - n, , 


So that for n = пу we have da(n) = 2 while don) > 2 for і < К. Since we can choose 
our primes anyway we like as long as p; # 2 this construction can be carried out in 
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infinitely many ways. If we ask for the smallest k for which daj(n) = 2 for all n then 
the answer is that such a k does not exist. 


Now to the problem. Factorizations and applications of (1) have been used to calculate 
k as a function of n for nx100, table 9. This does not tell us much about the genera] 
behaviour of dan). Table 10 provides some interesting cumulative statistics for n < 
10°. No more than 6 iterations are required for any n < 105. It seems strange that the 
column Юг К = 7 remains empty for n < 10? and n < 10? in particular in view of the 
regular behaviour in the column for К = 3 and our previous observation that К can be 
arbitrarily large for properly chosen sufficiently large n. This calls for further study. 


Table 9. К as a function of n for n < 100 


з 
x 
3 
x 
3 
x 
з 
< 
Ё 


1 
4 
1 
3 
3 
2 
1 
4 
] 
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Let p# denote the product of all prime numbers less than or equal to p and consider 
the largest number r of distinct prime numbers which are needed to construct any 
integer < 10° ie. p# < 10° < p.,&. With these primes consider all possible 
constructions 


293P.. p > 10° (2) 


This does not mean constructing all n < 10° but it does mean arriving at a structure 
into which all prime factorizations of n < 10° fits. This will be so because any number 
X 10° not produced by (2) will have fewer prime factors and smaller powers than one 
or more of the integers produced by (2). To illustrate this let's look at the case n < 
100. We have 2-3-5 < 100 < 2-3-5-7. We will therefore consider all possible 
construction 2*-3.5* < 100. These are obtained for a < 6, В < 4 and є < 2 resulting in 
table 11. 


Table 10. Number of iteration К required to arrive at ак(п)=2 







144 
1229 33 3444 3181 2108 4 
9592 79 34429 30466 24839 
271971 






о Ооо +N 


N 





1 
2 
3 
4 
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Апу number < 100 corresponds to one or more of these structures, for example 77 = 
7-11 corresponds to 1 1 0, 1 0 1 and O 1 1. This means that d(n) can only assume 
values listed in table 11 for n < 100. The above scheme has been computer 
implemented for n < 10”. The result is shown on table 12 and figure 1, which gives a 
clear picture of the overall behaviour of (п). 


Finally we will be able to say something about da(n). For n < 10° we have d(n) < 
6720<10*. From table 10 it is seen that К < 6 for < 10° and we therefore conclude 
that 


k € 7 forn < 107. 
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Table 12. Largest value of d(n) for п < TO К < 12 


k Largest Number of Number of Prime comb. Corresponding log{d) 
d d values combinations for largest d n 

1 4 4 4 11 6 1.3853 
2 12 11 34 121 90 2.4849 
3 32 22 141 3111 840 3.4657 
4 64 38 522 311110 9240 4.1589 
5 128 60 1848 3311010 98280 4.8520 
6 240 94 6179 4211101 942480 5.4806 
7 448 135 20198 63111100 8648640 6.1048 
8 768 190 42950 33211110 91891800 6.6438 
9 1344 266 133440 621111110 931170240 7.2034 
10 2304 359 399341 5312111100 9777287520 7.7424 
11 4032 481 783061 6322111100 97772875200 8.3020 
12 6720 626 2309712 64211111100 963761198400 8.8128 

107 

8T 

61 

щі 
44 
2+ 





Diagram 1. Largest value of In(diny) for n < 10. 


(b) Let о(п)= Уа and m a given positive integer. Find the smallest k for 
din.d>0 


which сої... 012)...)) 2m after k iterations, i.e. 904,2) 2m. 


Clearly К is a function of m. It is a stepwise increasing function. The first six 
iterations have been used to illustrate this for the interval 2 < m < 25 in diagram2. 


However, a far more interesting function to study is the inverse of k(m). This function 


m(k) is growing so rapidly that numerical results are difficult to interpret and 
represent. A better picture of the behavior of ск(2) is obtained by studying the function 


КК) = In(myk 
This function is represented in diagram 3 for the interval 1 < К < 100. After going 
through an interesting minimum for small values of k the curve flattens out. It seems 
to remain downwards convex. Does it approach a horizontal asymptote? 


Finally, a few iteration results (km): (1,3), (2.4), (3,7), (4,8), (5,15). (6,24), (7,60), 


(100,2972648508456959686477689735325484246606843303655482359755571200) 





10 15 20 25 


Un 


Diagram 2. k as a function of m. 
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Diagram 3. fík)-In(m//k. 
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1. F. Smarandache, Unsolved problem: 52, Only Problems, Not Solutions, ISBN 1-879585- 
00-6, Xiquan Publishing House, 1993 


58 


Chapter IV 


Diophantine Equations 


1. Some Thoughts on the Equation |y?-x*]-k 


In his book Only Problems, Not Solutions Smarandache formulated unsolved problem 
#20 as a conjecture rather than a problem: 


Let К be a non-zero integer. There are only a finite number of solutions in 
integers р, а, x, v, each greater than 1, to the equation x? - v! = К! 


In this study we will only consider the equation for pzq. By writing the equation in the 
form |х? - y?| = К we only need to consider cases where p > q. 


For a given set of parameters (р,4,К) it would then be desirable to list this finite 
number of solutions (x, y). However, if this were possible if would probably already 
have been done. So Smarandache's statement is likely to be based on statistical 
evidence rather than analytical reasoning. It is mainly from the statistical point of view 
we will study this equation. The parameters will be restricted to k<200 and р<9. As in 
most studies of this nature Ubasic provides the most effective computer language. All 
solutions where x«100 and у<100 can be churned out in a couple of seconds. To go 
any further we need a general approach to avoid running through meaningless search 
intervals. Consider 


x=4] tk 


For sufficiently large v only x = E гу" Тох = L4 Гу" J can produce solutions 
corresponding to 


1 Smarandache adds that “For k-1 this was conjectured by Cassels (1953) and proved by Tijdeman 
1976)." (Gamma 2/1986) 
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Diophantine Equations 
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ке ر‎ 1-и (1а) 
kev у (15) 


It is easy to imagine from (1a) and (1b) that the number of solutions thin out rapidly as 
we increase y. Let's illustrate this for p=3 and q-2 by looking at the number of 
squares s which fit into the interval between (n-1) and (п+1)?, ie. integers s for 
which 


or 


(n-1)¥<s?<(nt1)? (2) 


Let Smin be the smallest and Smax the largest s which satisfies (2). We can then 
calculate the ratio f between the number of squares between (n-1) and (n*1) and the 
difference between these cubes, i.e. 


f PE NE 
(полу «(поту 

То have а solution we must have an s such that |s? - n°] = К. The smaller f is the 
smaller is the chance for this to happen for our very limited range for К. Let sj? be the 
largest square which is smaller than n? and s; be the smallest square which is larger 
than n? then s; and s; give an indication of the behaviour of |s? - n)| . These two Ways 
of looking at the problem are displayed in table 1 for a Sequence of values of n which 
have been chosen so that n is neither square nor cube. 


Table 1. Frequency of squares 








19 4 3.95.102 
50 316 391 1.40.1033 











5-10? 14761 7600 4.40.1 0:5 
5-103 430916 276191 1.41۰106 
5-104 2515600 19845079 4.47-108 
5-105 287598881 419507900 1.41-10-9 
5-106 11203852544 11156827231 4.47.10 
5-107 513527672836 193579108351 1.41.1012 






5.108 1151976475001 21208703299996 4.47-10 4 





Table 2. Largest y 








С 

24 3 2 8158 736844 
199 4 2 10 99 
28 4 3 15 37 
60 5 2 76 50354 
24 5 3 4 10 
13 5 4 3 4 
127 6 2 4 63 
37 6 3 2 3 
104 6 4 3 5 
32 6 5 2 2 
95 7 2 6 529 
10 7 3 3 13 
95 7 4 6 23 
96 7 5 2 2 

64 7 6 2 

161 8 2 3 

8 3 2 

8 4 2 

8 5 2 

8 6 2 

8 7 2 

9 2 3 

9 3 2 

9 4 2 





Although the search for solutions was extended to у- 10000 no solutions were found 
for the following values of k: 


6, 14, 21, 29, 34, 42, 43, 50, 52, 58, 59, 62, 66, 69, 70, 75, 78, 82, 85, 86, 91, 102, 
110, 111, 114, 123, 125, 130, 133, 134, 146, 149, 150, 158, 160, 165, 173, 176, 177, 
178, 182, 187, 189, 195. 


No solutions were found for (p,q)=(9,5), (9,6), (9,7) and (9,8). 


Computer search: Calculations were carried out in Ubasic for y<10000, k<200 and 
р<9. Most solutions occur for small values of y for which we cannot use ће CEL ( 1) 
and the FLOOR (| J) functions. The ROUND function has been used instead. 
Calculations are carried out with real numbers to 19 decimal places (POINT 9 in the 
Ubasic language).To safegard against “near integer" solutions a “proposed” solution is 


61 


recalculated with integers in а subroutine. А simple version of the program is given 
below: 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
40 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 





point 9 

У1=10000 

for R%=3 іо 9 

for S%=2 to Ё%-1 

for K%=1 to 200 

Y=2 

while Y<Y | 

if (YAR%-K%)<1 then X1=(K%-YARB)A(1/S%):goto 100 
Х1=(У^В%-К%)^ ^(1/596) 

X-round(X1) 

if abs(X-X1)«10^(-30) then gosub 210 
endif 

Х1=(Ү^ЛЕ%+К%)^(1/$%) 

X-round(X1) 

if abs(X-X1)«10^(-30) then gosub 210 
endif 

inc Y 

wend 

next:next:next 

end 

if abs(YARKR-XAS%}<>K% then goto 240 
if X<=1 then goto 240 

print R%,S%,Y,X,K% 

retum 


The number of solutions for each set of parameters 15 given in table 3. The largest 
value of у which occurs in a solution for each parameter set (p,q) is given in table 2, 
which confirms the indications for the rarity of solutions for large y given in table 1. 


The largest number of solutions (11) occur for К=17. Several of these are due to the 
fact that no distinction is made between (x°)? and (x^. These solutions are displayed 
in table 4. A limited search (у<100) was carried out for 10<р<20, k«200. Only two 
solutions with yz2 were found. These results are shown in table 5. 


Conclusion: Smarandache's conjecture is well supported by the numerical results 
obtained in this study. The number of solutions diminish rapidly with increasing y, p 


and q. 
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Table 4. Solutions for kz17 
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2 
2 
2 
2 
2 
3 
2 
2 
2 
2 
2 
2 
3 
7 
2 
2 
2 
2 
3 
2 
2 
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Table За. The number of solutions to |ye-xa|-k 


p 344555666677777888888% 
ка 2 2 32 3 4 2 3 4 5 2 3 4 5 6 2 3 4 5 6 7 2 Sum 







Zu ca КУ 
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Table 3b. The number of solutions to |ye-xaj=k 


p 3445556666777778888889 
Ка 223 2342345 2 3 4 5 6 2 3 4 5 6 7 2 Sum 


2 
1 
6 
6 
4 
1 
3 
3 
5 
4 
3 
5 
2 
2 
64 ] 1 1 1 5 
4 
1 
4 
5 
2 
6 
2 
1 
1 
3 
3 
2 
2 
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Table Зс. The number of solutions to |УР-хч|=К 


р 344555666677777888888999 
ма 22 3 23 4 2 3 4 5 2 3 4 5 6 2 3 4 5 6 7 2 3 4 Sum 


Kiisa 
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3 
1 
1 
1 





Table За. The number of solutions to |ye-xs|-k 


р 3 4 4 5 5 5 6 6 6 6 7 7 7 7 7 8 8 8 8 8 8 9 9 9 
ka 2 2 3 2 3 4 2 3 4 5 2 3 4 5 6 2 3 4 5 6 7 2 3 4 Sum 


ca 


EAE ET КУ 
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2. The Equation 7(p*+q‘+r*+s*+t‘ - 5(р?+д?+т? +s+} + 90pqrst = 0 


This equation first appeared in the Numbers Count column of the Personal Computer 
World in March 1986. At that time the author found two until then unknown solutions 
(P.4,1,8,t)}=(87,42,6,3,3) and (p,q,r,s,t}=(99,97,39,13,2). This took approximately 75 
hrs on an 8086 processor running a 4.7 Mhz. Now, in 1996, these results were 
reproduced using the same program in 21 minutes on a Pentium 100 Mhz. The method 
used is reproduced below together with a number of new solutions. 


Consider the function 
F(p.q,r,s;t) = "(рани нн) - S(p^-q^ er^ es^ ey! + 90pqrst 
This function is invariant under exchange of variables. This makes it useful to study 


the function while all variables except one is kept constant. Denote this variable x. 
After some elementary algebraic manipulations we have 


F(...,x+*1,...)= F(...,x,...) + G(x) - 10C-H(x) + D/x 


where 
G(X) = 28x? + 22x? + 8x +2 
НО) = 2x + 1 
С = p^eq'«r^es^ et? with one of the constants replaced by x 
D = 90pqrst with the same constant as above replaced by x 


This permits us to calculate F(...,.x+1,...), G(x), H(x), C and D. When, in the next 
step, one of the unknowns p,q,r,s,t is increased by one the following replacements 
must be made 


C:i= C+ Hex) 
D = IXx+1 yx 


Without imposing any restrictions we can assume p» >г > >1. Fora given value 
of p the search will be conducted for descending values of the other unknowns. For 
given values of p, q, T, s consider the function 


p(t) = 7(ен) - 5(a + 2) +bt 


where e = p*tq‘tr*ts*, a = pqr +5? and b = 90pqrs 
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We have 
g(t) = 8? - 20at +b 
g (t) = 24€ - 20a 
g (t) = 0 for tm = (5a/6) which gives g°mn = b - 163 


Case 1: g’min > 0. Since g’(0)>0 it follows that g'(t)»0 for t»0. This means that g(t) is 
an increasing function for t>0 If we have found t; such that g(t;)»0 then g(t)>0 for all 
t^t; and the search can be interrupted for t-t;. 


Case 2: 2’minS0. For t>tm the function g(t) is convex and a value t=t; will be found for 
which the function is positive and increasing. The search is stopped for t>t;>tm. 


The above method has been used in a computer program written in Ubasic. 
Implementation on Pentium 100 Mhz computer has revealed a few new solutions. A 
complete list of all solutions for p<400 is given in table 4. 


Table 4. Solutions for px400 


# р а r 5 t 
1 1 1 1 1 1 
2 2 1 1 1 1 
3 2 2 ] 1 1 
4 3 3 2 ] 1 
5 4 2 1 1 1 
6 6 3 2 ] 1 
7 7 7 4 2 1 
8 17 7 7 1 1 
2 59 47 19 7 2 
10 87 42 6 3 3 
11 99 97 39 13 2 
12 124 63 42 17 ] 
13 127 47 34 2 1 
14 189 87 27 3 3 
15 286 154 11 11 11 
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3. The Equation у = 2-х1х2 ... xy +1 
Conjecture: 


Let k22 be a positive integer. The diophantine equation: у = 2-xix; ... Xx + 1 
has a infinitely many solutions in distinct primes y, X1,X2, ... Xk. 


This is unsolved problem number 11 in Smarandache's book Only Problems, Not 
Solutions. The word distinct has been added by the author. The purpose of this study 
is to see ‘how stable’ this conjecture is. This is done through a computer analysis of 
all possible solutions for y < 10°. (A very thin layer when surfing the ocean of integers 
but big enough to take a bit of time on the computer when it comes to calculations -in 
fact more than 100 hrs). From the computational point of view there is no reason to 
exclude k=1. The interval 0<у<10? is divided into 10 sub-intervals of equal length; 


Interval #1: 0<у<10® 
Interval #2: 10®&<y<2-10° 
Interval #3: 2-10%&<y<3-108 
Interval #10 9-10%&<y<10? 


The endpoints are excluded since these do not contribute to the number of solutions. 


Consider 
t=(v-1)/2=-x1x2 ... Xk а) 


The task is to identify sequentially all square free numbers пс 10?, For each square 
free number with k distinct prime factors we calculate the corresponding number y 
and test whether it 15 prime or not. The number of primes is denoted шк and the 
number of square free numbers 1s denoted ny. my and ny are recorded for each interval 
and the frequency of solutions Fy-my/ny is calculated. The result is shown in table 5. 


The same result are shown in diagram 1, which conveys a good visualization of a 
result obtained through surfing on a small area of the ocean of integers. 


Let's make a few observations: 


Why the irregularities for К-8 and k=9? The smallest square free integer with К prime 
factors is руй where px is the kth prime. For k=8 and k=9 this means there can be no 
solutions for y<2-19#+1=19399381 and y<2-23#+1=446185741 respectively. The 
samples for k=8 and k=9 are therefore too small to give a true picture - randomness 
takes precedence over mass behaviour. 
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Table 5. Frequency of solutions 






0.1282 
0.0701 0.0817 0.0967 0.1159 
0.0683 0.0794 0.0937 0.1118 
0.0672 0.0778 0.0919 0.1095 
0.0662 0.0769 0.0903 0.1080 
0.0655 0.0760 0.0895 0.1063 
0.0653 0.0752 0.0887 0.1054 
0.0646 0.0748 0.0878 0.1046 
0.0641 0.0744 0.0873 0.1036 
0.1030 


о MANA © М WKH м 





~ 
e 








каб 


0.2607 
0.2810 
0.2622 
0.2252 
0.2298 
0.2496 
0.2225 
0.2359 





0.0000 
1.0000 
1.0000 
0.5000 
0.6667 








Interv al 





Diagram 1 .Frequency of solutions 
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With the previous observation in mind let's compare the frequencies for the first and 
the last interval. Diagram 2 shows Ек for the intervals 0«y«10? and 9-105«y«10? . 


Frequency of solutions for intervals 1 and 10 























Diagram 2 


Observation: When not too close to p# the frequency smoothly increases with 
increasing k. This is a good support for Samrandache’s conjecture. 


Let’s now take a closer look at how the frequency of solutions behave for distinct 
values of k. Diagram 3 shows that the frequency decreases slightly for all k as larger 
integers as included. Is there an asymptote for each k?. The frequency of solutions 
increases as we increase k. It should be noted that it is the ratio between the number 
of solutions and the number of square free numbers in an interval determined by t in 
(1) which depicted. This is of course different from the number of square free 
numbers in the interval for у because when у runs through the interval a<y<b then 
we consider square free numbers s which obey 


(a-1y2«s«(b-1y2 
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As a bi-product to this study we have information on how many square free integers 


with 1, 2, ...., 9 prime factors there are in our different intervals. This is shown in 
table 6, where, as we have seen, we now only have five intervals. 

Interval £1: 0<s<108 

Interval #2: 105«s«2-10* 

Interval £3: 2.108<5<3.10% 

Interval #4: 3.105«s«4-105 

Interval £5: 4-105«s«5-108 


Frequency of solutions by intervals 


4 6 


Interv al й 





Diagram 3. 


Table 6. Number of square free numbers 






k=8 k=9 





k=7 














5761455 17426029 20710806 12364826 3884936 605939 
5317482 16565025 20539998 13029165 4476048 799963 63642 1409 
5173395 16270874 20457818 13243252 468954! 879765 76114 2060 
5084001 16085983 20402004 13374830 4825914 932513 84968 2560 
5019541 15951738 20359052 13468885 4926227 972398 91767 3005 





о с N O о 
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Diagram 4 shows a cumulative representation of the number of square free numbers 
with k71. 2, 3, ... 9 prime factors. A square free number with only one prime factor is 
indeed a prime number so the first column in table 6 shows the number of prime 
numbers in the interval. 


Cumulative statistics of square free numbers 


-4 k-5 k= k-7 k8 k-9 





Diagram 4. 

The illustration on the cover is another version of the above diagram for п<5-107. 
Finally, a problem often leads to other problems. Let's go back to iterations. In the 
introduction to chapter III it was said that iterations result in an invariant, a loop or 


divergence. Is there really no other possibility? Let's look at this: 


We define щ- = 2-uy*l where u; is a prime number. их. will have the 
property of being a prime number or not being a prime number. 


In this case we can give an explicit formula for ux-; in terms of u;. One easily finds 
ux; = 2(uy+1 1. 
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How do we characterize the result of iterations in relation to the property in which we 
are interested? Indefinite? 


How many times can we iterate щ- = 2-ux+] preserving primality? For ш=305192579 
it is eight times resulting in the following series of nine primes: 


305192579, 610385159, 1220770319, 2441540639, 4883081279, 
9766162559, 19532325119, 39064650239, 78129300479 


Which is the first series with 10 terms - or 11? Maybe we need some deep sea diving! 
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The Smarandache Сей Function 


Definition: For a positive integer n the Smarandache ceil function of order k is 
defined through ' 


$«(п) = m where m is the smallest- positive integer for which n divides m“. 


In the study of this function we will make frequent use of the ceil function defined as 
follows: 


| [х1 = the smallest integer not less than x. 


The following properties follow directly from the above definitions: 


1. Ѕ(п)=а 
2. SKp™ = реч! for any prime number р. 
3. For distinct primes р, q, ...r we have S(p^q*--0) = ا ...و ام‎ 


Henry Ibstedt 


$ 9.95 


